Large adiabatic geometric phase in rotating single atom-scale diamond defects 
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We demonstrate that the internal magnetic states of a single nitrogen- vacancy defect, within a 
rotating diamond crystal, acquire geometric phases. The geometric phase shift is manifest as a 
relative phase between components of a superposition of magnetic substates. We demonstrate that 
under reasonable experimental conditions a phase shift of up to four radians could be measured. 
Such a measurement of the accumulation of a geometric phase, due to macroscopic rotation, would 
be the first for a single atom-scale quantum system. 
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The quantum geometric phase is at the core of our 
understanding of the non-intuitive quantum view of the 
world, with historical origins going back to the original 
paper by Aharonov and Bohm [h, in 1959. However, 
while there have been many experimental demonstrations 
on ensembles of quantum systems, no measurement of a 
quantum phase on an individual quantum system, under- 
going macroscopic rotation, has been performed. In this 
paper we analyze the appearance of the geometric phase 
in defects in diamond and show how the measurement of 
the geometric phase of a single quantum spin, undergoing 
macroscopic rotation, is now in reach. 

In 1984, M.V. Berry developed an elegant and pow- 
erful mathematical framework which established the 
Aharonov-Bohm effect as just one instance of a far more 
general class of phenomena [2] ■ Berry considered the evo- 
lution of a system under a Hamiltonian which is adiabati- 
cally changed over time. He showed that the state of such 
a system acquires a phase which is geometrical in nature. 
The phase depends only on the system's path in param- 
eter space, specifically the flux of some field enclosed by 
that path. The field in question is a gauge field, akin 
to those found in quantum field theories, which arises 
naturally in Berry's formalism. 

Berry's work has since been applied to a diversity of 
phenomena, which can be broadly grouped under the 
umbrella of geometric phases or topological phases 0E]- 
Specific instances of these geometric phases include the 
following: various analogues of the Aharonov-Bohm ef- 
fect such as the Aharonov-Casher [5], Casella [5] and the 
He-McKellar [7] phases; the rotation of the polarisation 
of light in twisted optical fibres, which was recognised 
by Pancharatnam well before Berry's paper [5]; the so- 
called molecular Aharonov-Bohm effect which introduces 
a gauge field to nuclear degrees of freedom in the Born- 
Oppenheimer approximation to molecular dynamics [S]; 
and even the dynamics of classical systems such as low 
Reynolds number hydrodynamics |10j . 

Geometric phases have also proved a fruitful avenue 
of investigation for mathematical physicists due to their 
rich topological properties and their close connection 




Figure 1: (a) Energy level diagram of the NV center, (b) 
Geometry of NV center. Defining the magnetic field direction 
of the microwave pulses as the z-direction, z is the instan- 
taneous direction of the NV axis, defined with respect to the 
lab-frame, unprimed coordinate system, by and <j>. 

with gauge theories of quantum fields. Examples include 
mathematically formulating geometric phases in terms of 
the holonomy of line bundles [llj and directly using the 
geometric phase to help explain fractional statistics in the 
quantum Hall effect [T2] and the origin of Wess-Zumino 
terms in theories of quantum chromodynamics |13) . 

Despite this wealth of applications and observations, 
to date, only a few experimental observations of geomet- 
ric phases due to mechanical rotation have been made 
[HI ITS] . These measurements have been on ensembles 
of 35 C1 [T3] and 131 Xe [15] nuclear spins. However, no 
such measurement has been performed on an individual 
quantum state. 

The diamond nitrogen-vacancy (NV) system presents 
itself as an excellent tool for studying geometric phases. 
The electron spin is the canonical quantum system and 
the NV center offers a system in which a single spin can 
be initialised, coherently controlled, and measured. It is 
also possible to mechanically move the diamond crystal, 
and the NV with it, about some cyclical macroscopic 
trajectory within the spin coherence lifetime. As such 
this is an ideal system to investigate the accumulation 
of geometric phase, due to macroscopic rotation, on an 
individual quantum state. 
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The NV defect has a spin triplet ground state with a 
2.88 GHz zero-field splitting between the m = state and 
the degenerate m = ±1 inert states [Fig. 1(a)]. Optical 
excitation with 532nm light can pump the defect into 
the m = state and allows the population of the ground 
state to be read, since the m — state produces more 
fluorescence than the m — ±1 states [THl [T7]. The ef- 
fective Hamiltonian of the ground state, ignoring crystal 
asymmetries and hyperfine effects, is 



h = Id sl 



B S. 



(1) 



The first term is the zero-field splitting of the NV system 
itself, where D is the zero-field splitting strength. It is 
this term which makes the crystal's orientation crucial. 
It defines a quantisation direction z' (we reserve z for 
the lab-frame coordinate) which lies along the axis con- 
necting the nitrogen atom to its adjacent vacancy [Fig. 
1(b)]. The second term is the usual Zeeman splitting in- 
teraction with a magnetic field. The magnetic field is B, 
and g ~ 2 and \xb are the Lande g factor and the Bohr 
magneton respectively. 

To evaluate the geometric phase for the NV center the 
instantaneous eigenstates of the Hamiltonian in terms 
of the adiabatically varied parameters 9 and <f> [see Fig. 
1(b)] need to be determined. Writing the zero field split- 
ting Hamiltonian as 
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The geometric phase, given initialisation into \m) z >, 
after the crystal has been rotated along some trajectory 
P, defined by <f> and 6 is given by 



m (1 — cos 6*) d(j>. 



(4) 



For a closed loop in parameter space, this is simply the 
solid angle enclosed by the trajectory of zf , see Fig. 2(a). 
In the derivation of Eq. Q we have assumed that the 
evolution of the system is adiabatic: +siv? 9(j) 2 <C 2D 2 . 
This implies that the adiabatic approximation is valid 
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Figure 2: (a) Proposed Ramsey geometry. The crystal is at- 
tached to a rotating spindle with the NV (z') axis of one of its 
NV centers at an angle 9 to the spindle (z) axis. Rabi pulses 
are produced by a microwave field linearly polarised with its 
magnetic field oscillating along the z-axis. The magnitude of 
the geometric phase after a complete rotation of the spindle 
is given by the solid angle subtended by the z'-axis. (b) Pro- 
posed spin echo geometry. The spindle is at an angle #o to the 
fixed z axis of the microwave field, and the NV axis is perpen- 
dicular to the spindle axis. The actual path of the NV («') 
axis encloses a solid angle of 2ir, which is not observable. The 
7r-pulse in the spin echo geometry, however, rectifies the al- 
ternating Berry phase accumulation, producing a total phase 
of 4#o after a full spindle rotation, the solid angle enclosed by 
the effective trajectory shown. 

provided the angular velocity of the crystal is much less 
that the zero-field splitting frequency, which is true for 
the cases we consider below. 

There is a gauge degree of freedom in defining the in- 
stantaneous eigenstates \m) z i. The geometric phase for 
a closed loop trajectory is independent of the choice of 
gauge. However, we are interested in a geometric phase 
for trajectories which are not closed loops and must then 
be careful in choosing our gauge. 

The geometric phase is observed through interaction 
with a microwave magnetic field Br and it is the phase 
difference between the NV center and the microwave 
field, as seen by the NV center, which is the observ- 
able quantity. Consider a linearly polarized microwave 
field tuned to the zero-field splitting transition, oscil- 
lating along the fixed z axis, Br = Bncos(ujt)z. The 
Hamiltonian for the interaction of this field with the NV 
center is 
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where the matrices are expressed with respect to the 
\m) z ' basis as defined by Eq. ([3]). The approximation ne- 
glects the term proportional to S z > and is valid for weak 
microwave fields: g^isBji/h <C oj. 
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Equation. ([5]) depends on both the polar 9 and az- 
imuthal <f> angles. The dependence on the polar angle is 
simply a matter of effective strength, i.e. the effective 
microwave field strength experienced by the NV center is 
just BftSmO. The dependence on the azimuthal angle is 
more interesting: it has the character of a phase. The ef- 
fect of a microwave pulse on state \ip) when the NV axis is 
at some azimuthal angle <j> = 4>q is the same as the effect 
when = and the eigenstate is modified by a phase fac- 
tor exp(— im z i(j)o). Since the linearly polarised microwave 
field can be decomposed into two counter-rotating fields, 
it should be no surprise that the effective phase of the 
NV center should be so closely linked to the microwave 
field's angle relative the the NV axis. 

This phase-like dependence of the interaction Hamilto- 
nian on the NV orientation can be eliminated via a gauge 
transformation, of the \m) z i basis, of the form 

\m) z ,^e^ 6 ^\m) z ,. (6) 

Choosing f m = —mcj) eliminates the ^-dependence of the 
interaction Hamiltonian [Eq. §5§) as desired. Using this 
gauge transformed eigenstate basis the geometric phase 
becomes 

$ = J m cos 9d(j). (7) 

The final state produced by the microwave pulse now 
depends only on the explicit phase given by Eq. Q. 

Because the diamond NV electron spin has such long 
coherence times it should be possible to observe a geomet- 
ric phase by mechanically spinning a diamond crystal. 
The set-up we propose is very similar to the Aharonov- 
Casher measurement proposed in Ref. [18]. A diamond 
crystal is mounted on a spinning spindle. Optical ini- 
tialisation, coherent microwave control, and fluorescence 
detection are used to measure the phase evolved after a 
certain rotation angle. Below we consider two possible 
scenarios for measuring the geometric phase: a relatively 
simple Ramsey geometry [Fig. 2(a)] and a spin echo pulse 
sequence [Fig. 2(b)], enabling a significant increase in the 
maximum geometric phase observed at the cost of intro- 
ducing additional n microwave pulses. In each case it is 
assumed that the pulse microwave pulse durations and 
readout times are of order 50ns and hence the diamond 
can be considered to be stationary. Finally the relative 
sensitivity of the measurements is investigated. 

The Ramsey geometry is shown in Fig. 2(a), with the 
microwave field linearly polarised with its magnetic field 
pointing along the spindle (z) axis. The crystal itself is 
mounted such that the NV axis makes an angle 9 to the 
spindle axis. We assume that some mechanism lifts the 
degeneracy of the |±l)z' states. This could be mechanical 
strain in the crystal, a static magnetic field along z (weak 
enough not to mix the NV eigenstates) or a magnetic 
field rotating with the spindle, produced, for example, 
by a permanent magnet mounted on the spindle. 

For the Ramsey geometry as the azimuthal angle <f> of 
the NV axis passes zero, a ir/2 Rabi pulse is applied, 



tuned to the |0) z / — > \l) z > transition. As the spindle 
continues to rotate, a relative geometric phase evolves 
between the |0) z < and \l) z / states, given by $ = 4>cos9. 
After the spindle has rotated through some angle cf>o, a 
second tt/2 pulse is applied, converting the phase into 
a population difference, which is measured by 532 nm 
illumination and a fluorescence recording. A spindle ro- 
tation speed of / = 2kHz and a lOfis inhomogeneous 
broadening time (T% ) [TS] would allow up to 20mrad of 
geometric phase to be measured. 

To extend the coherence lifetime of the NV electron 
spin a spin echo pulse sequence could be employed, hence 
enabling a larger geometric phase to be measured. To do 
this a different geometry is required, shown in Fig. 2(b), 
in which the spindle axis is placed at an angle 9q to the 
microwave field (z) axis, and the NV (z 1 ) axis is per- 
pendicular to the spindle axis. The resulting geometric 
phase, alternates, and for a complete rotation is zero. 

A spin echo control sequence, with tt pulses applied 
whenever the NV (z 1 ) axis is perpendicular to z, would 
rectify the alternating geometric phase, producing a total 
phase of $ = 4n9, where n is the number of complete 
spindle rotations. By extending the coherence time from 
T| to T2 « 2ms [20] a geometric phase of 4 radians could 
be produced. 

The measurements proposed above would consist of 
many repetitions of the pulse sequence to get an average 
signal. The sensitivity of the measurement of the geo- 
metric phase is determined by the Poissonian statistics of 
spin projection, photon emission, and photon collection. 
The uncertainty A$ in the measurement of a geometric 
phase is related to the uncertainty AS of the normalised 
fluorescence signal S, 



A4> = AS (-1 = 2AS. (8) 

The second equality arises because the normalised signal, 
a sinusoidal function of $, has a maximum gradient of 
1/2. By appropriately retarding the phase of the final 
it/2 pulse it can be ensured that the sinusoid is at its 
steepest point at the time of measurement. 

The normalised signal S is the number of photons col- 
lected over N r runs, normalised so that (S) =1/2 when 
the populations of |0) and |1) are equal. If each mea- 
surement of |0) or |1) corresponded to the emission and 
detection of exactly one or zero photons respectively, 
then the variance of S would be (AS) 2 = l/(2N r ). A 
more careful analysis, taking into account the statistics of 
(spontaneous) photon emission, imperfect detection and 
the nonzero fluorescence of the |1) state, modifies the 
variance of the normalised signal by a factor C 2 , giving 
(AS) 2 = l/(2C 2 N r ). The physical basis for the factor C 
(C R3 0.15 for typical experiments and C = 1 in the ideal 
case) is described by Taylor et al. [2"T] . 

The relative sensitivity for a series of geometric phase 
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measurements, using a single NV center, is 
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where Tm is the time to take a single measurement, Tt is 
the total averaging time. Expressing Tm in terms of the 
relevant de-coherence time (Tm — clT^*' ', where a > 1) 
the relative sensitivity becomes 
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Based on the proposed experimental parameters given 
above (a = 2), this corresponds to a relative uncertainty 
in the measurement of $ with (without) spin echo pulses 
sequences of 0.15 Hz" 1 / 2 (2 Hz" 1 / 2 ), or a 0.15% (2%) 
uncertainty after three hours. 

We have demonstrated that a geometric phase shift 
manifests between the internal magnetic states of a single 



nitrogen-vacancy defect, within a rotating diamond crys- 
tal. The measurement of such a geometric phase shift 
in a macroscopically rotating single atom-scale quantum 
object would provide a unique test of our fundamental 
understanding of quantum mechanics. As such we have 
demonstrated that the measurement of geometric phase 
shifts of > 1 radian in such systems is possible. The 
analysis presented above is not only important in terms 
of demonstrating geometric phase shifts in macroscopi- 
cally rotating quantum systems it also provides the basis 
for quantifying geometric phase shift effects in the use of 
nano-diamonds as high precision translational and rota- 
tional sensors 1221. 
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